Nonlinear Excitations, Stability Inversions and Dissipative Dynamics in 
Quasi-one-dimensional Polariton Condensates 
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We consider the existence, stability and dynamics of the ground state and nonlinear excitations, in 
the form of dark solitons, for a quasi-one-dimensional polariton condensate in the presence of pump- 
ing and nonlinear damping. We find a series of remarkable features that can be directly contrasted 
to the case of the typically energy-conserving ultracold alkali-atom Bose-Einstein condensates. For 
some sizeable parameter ranges, the nodeless ("ground") state becomes unstable towards the forma- 
tion of stable nonlinear single or multi dark-soliton excitations. It is also observed that for suitable 
parametric choices, the instability of single dark solitons can nucleate multi-dark-soliton states. 
Also, for other parametric regions, stable asymmetric sawtooth-like solutions exist. Finally, we con- 
sider the dragging of a defect through the condensate and the interference of two initially separated 
condensates, both of which are capable of nucleating dark multi-soliton dynamical states. 



I. INTRODUCTION 

An important recent development that has spurted a 
new direction for the physics of Bosc-Einstcin condensa- 
tion has been the observation of condensation phenom- 
ena for polaritons in semiconductor microcavities [l[ at 
much higher temperatures than ultracold atomic Bosc- 
Einstein condensates (BECs) [1, H|. In the setting of 
exciton-polariton BECs, the condensing "entities" are 
excitons, i.e., bound electron- hole particles. When con- 
fined, these develop strong coupling with light, forming 
exciton-photon mixed quasi-particles known as polari- 
tons 0| . The finite temperature formation of polariton 
condensates leads the quasi-particlcs to also possess a fi- 
nite lifetime: in fact, they can only exist for a few picosec- 
onds in the cavity before they decay into photons. Hence, 
thermal equilibrium can never be achieved and the sys- 
tem produces a genuinely far-from-equilibrium conden- 
sate, in which external pumping from a reservoir of ex- 
citons counter balances the loss of polaritons due to the 
above decay mechanism. Nevertheless, numerous key fea- 
tures of the supcrfluid character of the exciton-polariton 
condensates have been established, including the flow 
without scattering (analog of the flow without friction) 
[j| , the existence of vortices @ (see also Ref. 01 for vor- 
tex dipole dynamics) , the collective dynamics [8| , as well 
as remarkable applications such as spin switches Q and 
light emitting diodes floj operating even near room tem- 
peratures. 

The pumping and damping mechanisms associated 
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with polaritons enable the formulation of different types 
of models. One of these, proposed in Refs. [ll|, [l2[ 
suggests the use of a single partial differential equation 
(PDE) for the polariton condensate incorporating the 
above mentioned loss-gain mechanisms. This model fea- 
tures a localized (within a pumping region) gain and 
a nonlinear saturating loss of polaritons; these are the 
fundamental differences of this setting from the stan- 
dard PDE mean-field model, namely the Gross-Pitaevskii 
equation (GPE) used in the physics of atomic BECs 
H, EH- In another class of models, which has been 
proposed in Refs. [HHHI, the polaritons are coupled to 
the evolution of the exciton population; such models also 
display nonlinear diffusive spatial dynamics for the exci- 
tons. 

In this work, our aim is to consider the quasi-one- 
dimensional (ID) dynamics of polariton BECs and to 
illustrate their fundamental differences from the more 
standard alkali-atom condensates. At this point we 
should mention that polariton condensates considered 
so far (also experimentally) have been intrinsically two- 
dimensional (2D) [l7|. Nevertheless, using tight con- 
finement along one, transversal, direction (i.e., highly 
anisotropic variants of the traditional parabolic traps) 
we envision rendering the polariton condensate effectively 
ID [l8j . In addition, a broader perspective for our con- 
siderations is that understanding the nonlinear dynam- 
ics and pertinent phenomenology in the ID setting, may 
pave the way towards subsequently generalizing relevant 
considerations to the more realistic 2D case. The key 
phenomena that are reported herein are: a wide paramet- 
ric interval of destabilization of the fundamental nodeless 
state of polariton BECs; a partial stability within this in- 
terval of excited states, in the form of dark solitons; the 
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spontaneous production of higher excited (multi-soliton) 
states from lower ones or even from nodeless states; the 
production of dark soliton trains by dragging of a defect 
through the polariton BEC (cf. the recent relevant 2D ex- 
perimental results in Ref. [191]); and finally the formation 
of long-lived dark multi-soliton dynamical states through 
the interference of two separated polaritonic condensates 
in analogy with the atomic BEC case of Ref. [20| • 

The paper is organized as follows. In Section II we 
describe our model and setup, providing also a brief de- 
scription of our methods. Section III is devoted to our 
detailed numerical investigations, and Section IV con- 
cludes our work, including some suggestions for possible 
future studies. 



II. MODEL SETUP 

In our analysis below, we consider the modified com- 
plex Gross-Pitaevskii model developed in Refs. [TJ E3 
suitably reduced to one spatial dimension: 

iW = {-d 2 x + x 2 + |Vf + i [(*(*) - a|V| 2 )] } V, (1) 

where ip denotes the polariton wavefunction trapped in- 
side a ID harmonic potential, x 2 (the transverse direc- 
tion, perpendicular to x, corresponds to the tight trap- 
ping axis mentioned above). The differences of Eq. ([1]) 
from the standard GPE appearing in the physics of 
atomic BECs can be traced to the presence of (i) the 
spatially dependent gain term with 



x(x) = aQ(x m - \x\), 



(2) 



where is the step function generating a symmetric 
spot of "radius" x m and strength a for the gain and (ii) 
the nonlinear saturation loss term of strength a. Esti- 
mates of the relevant physical time and space scales, as 
well as physically relevant parameter values, are given in 
Ref. We should also note in passing that although 
our results below are given in the context of Eq. (JT|), wc 
have ensured that similar phenomenology arises in the 
model of Refs. [I3 - [T6j . for suitable parametric choices. 
In other words, the phenomenology that is reported in 
this work is generically relevant to (ID) polariton BECs 
independently of model specifics. 

In what follows, we will consider the stationary so- 
lutions of the quasi-lD model at hand, in the form 
ij)(x,t) = ipo(x) cxp(-ifit) where fi is the dimensionlcss 
chemical potential, and the stationary state ipo( x ) is gov- 
erned by the following ordinary differential equation: 

M) = j-^a + x 2 + |^o| 2 + i [(x{x) - frivol 2 )] | V'o- 

(3) 

Importantly, the additional condition 



needs to be enforced as a population balance constraint. 
This self-consistently selects the particular value of the 
chemical potential once the other parameters (i.e., a, a, 
x m ) are fixed. This is why some of our graphs of the solu- 
tion branches below will feature \x as a function of other 
solution parameters, such as x m . We note in passing the 
significant difference of this trait from the Hamiltonian 
atomic BEC case, where there exist monoparmctric fam- 
ilies of solutions as a function of /i. Once stationary so- 
lutions of the differential- algebraic system of Eqs. ©-(lU 
arc identified, their linear stability is considered by means 
of a Bogolyubov-de Gennes analysis. Namely, small per- 
turbations [of order 0(<5), with < 8 <C 1] are introduced 
in the form 

i>(x, t) = e~ iflt U (x) + 8(a(x)e lut + b* (x)e~ 1 ^ **) 

and the ensuing linearized equation are then solved to 
0(6), leading to the following eigenvalue problem: 



a(x) 
b(x) 



L\ L2 
-1*2 —L* 



a(x) 
b(x) 



(5) 



for the eigenvalue w and associated eigenvector 
(a(x) , b(x)) T , where L\ and L 2 are the following oper- 
ators: 



Li 



-M - + x + 2 ( X _ i(J )\i'o\ + ix(x), 



L 2 = (1 - icr)V> 2 . 

Once the stationary solutions are found to be linearly 
unstable (i.e., Im{w} 7^ 0), then the dynamical mani- 
festation of the corresponding instabilities is monitored 
through direct numerical simulations of Eq. JT]). In addi- 
tion, in what follows, we have considered dynamical sce- 
narios under which nonlinear excitations, such as single- 
or multiple-dark-solitons [l3| can arise in the context of 
polariton BECs. Such excited states have been exten- 
sively studied in the context of atomic BECs [2l| , while 
they have been amply considered in recent experimen- 
tal investigations in this context [2(| HH, HH- Addi- 
tionally, motivated by relevant studies in atomic BECs 
21J25L and recent experiments in polariton condensates 
5l Il9l. |26| , we consider the nucleation of dark solitons by 
a moving defect, modeled by a (localized) potential; the 
latter, is assumed to be produced by a narrow laser beam 
of Gaussian shape, namely: 



Vde 



y exp(-(x-i>t) 2 /e 2 ), 



(6) 



da; (x(x) 



a|^| 2 )|^ | 2 = 



(4) 



where Vo, v and e represent the amplitude, speed, and 
width of the potential, respectively. We consider the fixed 
point solution of the modified GPE (Q]) in the presence of 
this second defect potential (in addition to the harmonic 
trap x 2 ) at the center of the trap. Then we evolve the 
system in time starting from this solution, dragging the 
defect through the system. This is similar in spirit to the 
recent experiments of Ref. IS | (and to the recent the- 
oretical investigation of Ref. 27 ) . Finally, we study an 
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FIG. 1: (Color online) Spatial profiles of the densities, lipf, 
of nodeless states and single dark solitons (insets). In the top 
(bottom) panel, the parameter a (x m ) is kept fixed, taking the 
value a — 2 (x m = 3), while a = 3.5 in both cases. As seen 
in the top panel, the profiles do not change appreciably above 
a critical value of x m ~ 3: in fact, the profiles corresponding 
to Xm — 3 [dotted (black) line] and x m = 5 [dotted-dashed 
(green) line] are almost indistinguishable from each other. All 
quantities in this and in all subsequent figures are dimension- 
less. 



alternative proposal for dark soliton nucleation, akin to 
the interference experiments conducted for atomic BECs 
in Ref. j2fj| . whereby a central potential barrier separat- 
ing two polariton clouds is lifted allowing the two clouds 
to interfere; this process leads to the production of per- 
sistent dark solitons as well. 



III. NUMERICAL RESULTS 

We hereby explore the existence, stability and dynam- 
ical properties of the nodeless cloud (NC) as well as of 
excited states exhibiting a single-node, namely dark soli- 
tons (DSs). The above mentioned states (NC and DS) 
are the most fundamental nonlinear states of the system, 
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FIG. 2: (Color online) Stability domains of nodeless clouds 
(NC) and dark solitons (DS) for a = 0.35 and a < 3.5. Dark 
solitons are stable only in the area indicated by the arrow. 
Notice that there is a superposition between part of the dark 
soliton and nodeless cloud stability ranges. The (green) circles 
indicate the parameter locations for the spectra depicted in 
Fig. El 



whose profiles —for different parameter combinations — 
are displayed in Fig. [TJ In particular, the top panel of 
Fig. [U depicts the NC and DS profiles for a constant sat- 
uration coefficient a = 3.5 and constant pumping spot 
strength a = 2, and a varying radius of the pumping 
spot x m . We observe that, for given values of a and a, 
there is a critical value of x m above which the shapes 
of the NCs and DSs remain unchanged. In the exam- 
ple depicted in the top panel of Fig. [IJ the profiles for 
x m <^ 3 are indistinguishable from each other (see pro- 
files for x m = 3 and x m = 5). This effect is due to the 
fact that for large x m , the pumping spot covers a larger 
portion than the saturated size of the cloud when loss 
and gain are balanced. This cloud size is the analogue of 
the Thomas-Fermi radius for an atomic condensate. The 
saturation of the cloud size is equivalent to the saturation 
of the chemical potential /i as x m is increased, as shown 
below (see Figs. 0] and [5]). In the bottom panel of Fig. Q] 
we depict the NC and DS profiles for a constant pumping 
spot radius x m = 3 and a varying pumping spot strength. 
Note that, in this case, the size of the cloud continuously 
expands (in amplitude and width) with increasing spot 
strength. 

We now proceed to provide a characterization of the 
existence and stability properties of the NC and DS pro- 
files with respect to the various parameters at hand. In 
what follows, in order to offer a picture of the relevant 
parameter space, we have varied the gain parameters a 
and x m whereas the coefficient of the saturating nonlin- 
ear loss a = 0.35 has been kept fixed. Figure [2] depicts 
the existence and stability domains in the (a,x m ) pa- 
rameter plane with a = 0.35 fixed for both the nodeless 
cloud and the single DS that can be found as (numerically 
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FIG. 3: (Color online) Spectral plane for nodeless clouds (top 
row in each series of panels) and dark solitons (bottom row 
in each series of panels) with a — 0.35. The top series of 
panels correspond to fixed a = 1 and increasing values of 
x m as labeled whereas the bottom series of panels correspond 
to fixed x m = 3 and increasing values of a as indicated. The 
different cases correspond to the parameter locations depicted 
by the (green) circles in Fig. 
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FIG. 4: (Color online) Top panel: Chemical potential as a 
function of x m for a nodeless cloud with a — 1.0 and a — 0.35. 
The insets show the time evolution of the cloud for various 
values of x m . 
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FIG. 5: (Color online) Top panel: Chemical potential as a 
function of x m for a dark soliton with a = 1.0 and a = 0.35. 
The insets show the time evolution of the cloud for various 
values of Xm • 



exact up to a prescribed accuracy) fixed point solutions 
of Eq. ([3]). The NC and DS configurations exist for all 
parameter combinations, as it is the case for atomic con- 
densates. 

Nevertheless, as far as the stability and dynamical 
properties of NC and DS states are concerned, we can 
observe fundamental differences between the pumpcd- 
dampcd polariton BECs and atomic BECs. In partic- 
ular, the nodeless cloud (which was always stable in the 
Hamiltonian case of atomic BECs [HI) is now stable only 
below a critical value of the pumping spot size x m . On 
the other hand, also remarkably, even the single DS is 
stable only in a limited range, while it was always stable 
in quasi-ld harmonically trapped atomic BECs (see, e.g., 
Ref. [2(| and references therein). Moreover, in a very un- 
usual manifestation of stability inversions, not only can 



the nodeless state be stable while the dark soliton is not, 
but also vice versa: the state with a node can be stable 
while the one without a node is not. In Fig.[3]we show the 
details of the Bogolyubov spectra of both states. These 
showcase the dissipative nature of the dynamics being as- 
sociated with frequencies chiefly with negative imaginary 
part; moreover, they also illustrate the potential instabil- 
ities arising in the system through Hopf bifurcations and 
oscillatory instabilities associated with complex eigenfrc- 
quencies, or through zero crossings (and purely imaginary 
eigenfrequencies) . This second scenario only appears for 
dark solitons. From an intuitive viewpoint, this phe- 
nomenology can partially be understood on the following 
grounds. The condensate in the absence of the external 
driving has an intrinsic length scale selected by the trap 
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(and the chemical potential). The presence of the exter- 
nal forcing over the radius x m introduces an additional 
length scale competing with the former one. Hence, when 
this forcing becomes fairly (spatially) extended, it favors 
a spatially wider state. This is manifested through the 
instability of a group of "background" spectral modes 
(which arc not the lowest modes of the condensate close 
to the spectral plane origin) in the panels of Fig. [3] 
On the contrary, the dark soliton or multi-soliton states 
may become unstable through the same mechanism, or 
they may also become unstable through their "internal 
modes" [2(| which lead to the isolated instability through 
the zero crossing. 

In Fig. 2] we show a continuation of the nodeless state, 
for a fixed value of a = 1.0. In particular, we depict 
the chemical potential fi as a function of pumping spot 
The insets show the dynamical evolution of the 
fixed point solution for several illustrative values of x m , 
which reflect the different possible dynamical phenom- 
ena. The values shown are in the unstable domain of the 
NC. There are two types of behavior. The first behav- 
ior, for 3.1 < x m < 3.8, corresponds to the NC decaying 
into multiple-dark-solitons, by means of oscillatory tran- 
sients (even in the region of stability of the single DS). 
The spontaneous emergence of these states from a node- 
less one is a feature particular to polariton BECs, having 
no analog in the atomic BEC case. On the other hand, 
a second behavior, corresponds to values x m > 3.8: in 
this case, even though transient multi-soliton states still 
appear, they finally give rise to nearly "turbulent" non- 
linear dynamics of a "sea" of multiple DS states, which 
may (or may not, depending on x m ) settle on an asymp- 
totic multi-soliton state. 

On the other hand, we have also investigated the dy- 
namics of the fundamental (single) dark soliton in Fig. [5] 
We have found that, in their instability region, DSs decay 
towards the nodeless state, as expected, if the latter is 
stable. However, when the nodeless state is unstable, the 
dynamics is as follows: after a transient stage, a breath- 
ing multiplc-DS structure is formed, which may consist 
of 3, 4 or 5 DSs. Examples of such evolutions are shown 
in the insets of Fig. [5] In this case, the DS is stable in 
the range 2.1 < x m < 3.1. 

The effects of varying the parameters a and x m , are de- 
picted in Figs. [5] and [71 for fixed pumping spot strengths 
a = 2 and a = 3, respectively. A new dynamical fea- 
ture that arises for a = 3 and sufficiently large a and 
x m (see Fig. [JJ is the emergence of a highly asymmetric 
"sawtooth" structure (notice that for this value of a, no 
stable DSs are found). Figure [5] shows how these saw- 
tooth structures emerge from the NC when the latter is 
unstable. The top panel shows the the emergence of a 
breathing sawtooth structure from an unstable NC. No- 
tice that the asymmetry of this configuration is produced 
by the amplification of the asymmetric (small) perturba- 
tion added to the NC in order to manifest its dynamical 
instability. We have checked that this breathing behav- 
ior persists for very long times (larger than several thou- 
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FIG. 6: (Color online) Top panel: Profiles for the nodeless 
clouds and dark solitons (inset) for different values of a and 
fixed q = 2 and x m = 2.5. Bottom panel: Stability domain 
for nodeless clouds (NC) and dark solitons (DS) for a = 2. 



sands) without any apparent decay. This is due to the 
fact that, for these parameter values, the steady state 
sawtooth configuration is unstable (see the spectrum de- 
picted in the inset of the top panel in Fig. [8]) and thus 
the cloud cannot decay to this state. For other param- 
eter values, the emerging sawtooth structures settle to 
stationary configurations without any breathing as can 
be seen seen in the bottom panel of Fig. [H This case 
corresponds to parameter values inside the elongated is- 
land depicted in the bottom panel of Fig. [JJ where these 
sawtooth configurations are stable. 

The top panel of Fig. [JJ depicts examples of station- 
ary sawtooth structures. Similar to the saturation of the 
nodeless cloud size for large pumping spot size x m , we 
have also observed a saturation of the size for the saw- 
tooth structures for large x m (results not shown here). 
It is extremely interesting that the polariton condensate 
is able to support stable asymmetric sawtooth-patterned 
states. Such solutions are not possible in atomic BECs. 
Furthermore, it is important to mention that, because of 
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FIG. 7: (Color online) Top panel: Profiles of highly asymmet- 
ric "sawtooth" structures for different values of a and fixed 
a — 3 and x m — 5. Bottom panel: Stability domain for 
sawtooth structures with a = 3. 



the left-to-right symmetry of our system, the sawtooth 
states appear in pairs (left- and right-handed sawtooth 
structures). In this work, for simplicity, we chose to only 
depict one family as the other one is symmetric with re- 
spect to the center of the cloud. 

Finally, although our results already suggest that non- 
linear excitations in the form of dark solitons should 
spontaneously emerge in polariton BECs, we offer some 
alternative dynamical schemes for producing such exci- 
tations inspired by experimental realizations within their 
atomic BEC counterparts [2(| [25[ , which also appear to 
be within reach for the case of polaritons; see e.g. the very 
recent work of Ref . (l9j and references therein. One possi- 
ble nucleating mechanism for dark solitons is by dragging 
an obstacle — in the form of the potential of Eq. flip— suf- 
ficiently fast through the condensate (see Refs. [25[ and 
[l9| for relevant experimental observations in atomic and 
polariton condensates, respectively). Examples of this 
effect are shown in Fig. [5] For the relatively strong har- 
monic confinement considered in Eq. ((T|), we have found 
that at most two DSs can survive, due to their relatively 



FIG. 8: (Color online) Density plots showing the evolution 
of unstable nodeless clouds to a breathing sawtooth structure 
for g = 0.5 (top) and a non-breathing sawtooth structure for 
a = I (bottom) with x m = 5. The inset in the top panel shows 
the linearization spectrum of the unstable sawtooth steady 
state about which the system oscillates about for long times. 
The inset in the bottom panel corresponds to the linearization 
spectrum for the stable sawtooth steady state that the system 
asymptotes to. 



large size. In this case, as shown in the top panel of 
Fig. [9l although three (or even four) DS can be seen 
being emitted out of the defect path, through collisions 
they eventually decay down to two, that continue to in- 
teract. The DS nature of these structures can be seen 
through the phase jump of 7r shown (together with its 
profile) in the inset. Chains of dark solitons (alias "dark 
soliton trains" ) can be produced by this dragging defect 
mechanism if one chooses a weaker harmonic trap Q 2 x 2 /2 
[instead of x 2 as in Eq. (JTJ], with a trap strength il suf- 
ficiently small. If fact, as depicted in the bottom panel 
of Fig. IH1 a weaker trapping with O = 0.04 (correspond- 
ing to a considerably wider condensate) allows for the 
formation of a train of DSs that propagates initially in 
the opposite direction of the dragging defect. However, 
contrary to the Hamiltonian case of atomic BECs where 
the distance between the generated DSs appears to be 
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FIG. 9: (Color online) Nucleation of dark solitons by a drag- 
ging defect traversing a nodeless cloud. Depicted are the den- 
sity plots of the nodeless cloud as it is traversed by the drag- 
ging defect. The dashed line represents the trajectory of the 
defect. The case depicted in the top panel corresponds to a 
tighter confinement with trap frequency = 1, and parame- 
ter values a = 1.0, a — 0.35 and x m — 2.8 whereas Vo = 8, 
v — 0.9 and e 2 = 4. The bottom panel corresponds to a 
weaker harmonic trapping with trap frequency f2 = 0.04, and 
a = 0.1, a = 0.08 and x m = 44.5, while Vo = 1, v = 0.5 and 
e 2 = 0.16. The insets show the profile ([green] dashed line) 
and its corresponding phase ([red] solid line) at the indicated 
times. 



approximately constant [28j . here DSs interact strongly 
with each other and with the background cloud, setting 
it in oscillation. This might also be due to the unstable 
nature of the underlying NC cloud for this parameter set. 
Again the DS nature is visible from the tt phase jump at 
each of the nonlinear excitations shown in the inset. 

Another possible nucleation mechanism for DSs is by 
nonlinear interference (see, e.g., Ref. [2l[ for a discus- 
sion in the context of atomic BECs). This mechanism 
can be realized by splitting the condensate into two frag- 
ments by adiabatically inducing a potential barrier in the 
central portion of a stationary nodeless cloud, and sub- 
sequently releasing the fragments by suddenly removing 
the barrier that separates them. This "nonlinear inter- 
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FIG. 10: Nucleation of dark soliton pairs by nonlinear inter- 
ference of colliding condensate fragments. The plots depict 
the density of the evolution of a cloud created in a double 
well and allowed to interfere by removing the central barrier. 
Parameter values are: (top) a — 1.0, a — 0.35 and x m = 2.8 
whereas Vo = 17 and e 2 = 1, fi = 1; (bottom) a — 0.1, 
a — 0.08 and x m = 44.5 whereas Vo = 4 and e 2 = 16, while 
the trap frequency is O = 0.04. The insets show the profile 
([green] dashed line) and its corresponding phase ([red] solid 
line) at the indicated times. 



ference method" was used in atomic BEC experiments 
to demonstrate the generation of vortex structures (29j 
(see also theoretical work in Ref. [3(|) and dark solitons 
[20l [23l . HJ . Figure [TU] depicts two samples of this nucle- 
ation mechanism where two polariton condensates, ini- 
tially placed in a double well potential, are released in the 
harmonic trap and allowed to interfere. Again we show 
the scenarios of strong (top) and weak(cr) harmonic con- 
finement (bottom). We considered an initial situation 
where a barrier of the same shape as (but wider than) 
the obstacle considered before [cf. Eq. [5])] is superposed 
at the center of the harmonic trap. At t = the bar- 
rier is removed, thus allowing the two portions of the 
condensate to mix and interfere. As can be seen in the 
top panel of Fig. [TUl a pair of dark pulses is formed; but 
contrary to the case of the soliton trains nucleated by 
the dragging defect, we observe that the DSs nucleated 
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by the nonlinear interference mechanism move relatively 
fast away from each other towards the edge of the con- 
densate. Notice that, as again shown in the inset, the 
phase jumps of approximately n indicate that these are 
pairs of genuine DSs. The bottom panel shows the case 
of weaker confinement where in this case more than one 
pair of nonlinear excitations is produced by the nonlinear 
interference. 



IV. CONCLUSIONS 

In this work, we considered and studied a com- 
plex Gross-Pitaevskii model describing the quasi-one- 
dimensional dynamics of polariton condensates. Our mo- 
tivation was to understand the fundamental differences 
between polariton condensates and their atomic (Hamil- 
tonian BECs) counterparts. We found that the specific 
pumping and damping terms that have been argued as 
being relevant to polariton condensates case offer a wide 
range of unexpected features when compared to atomic 
condensates. 

The fundamental nodclcss state (regarded as the 
"ground state") of the system was found to become un- 
stable through a variety of mechanisms, while excited 
states — in the form of single- or multiple-dark solitons — 
were found to result from the instability of the node- 
less state. The fundamental excited state, namely the 
single dark soliton was also found to be subject to in- 
stabilities, leading to either spontaneous formation of 
multi-dark soliton states or even to emergence of a "dark- 
soliton-turbulence" (when highly unstable). All these are 
dynamical manifestations which significantly distinguish 
the polaritonic case from the atomic BEC variant of the 
problem, yet we have intuitively attributed them to the 
emerging competition of length scales (among the intrin- 
sic length scale of the trapped system and the length scale 
of the applied forcing). We also observed the emergence 
of stable asymmetric sawtooth-like configurations which 
are not present in atomic BECs. 

Finally, other techniques, namely nucleation of dark 



solitons by dragging an obstacle through the condensate 
and through nonlinear interference of colliding conden- 
sate fragments, have been studied. These were adapted 
from the atomic condensate case in order to produce fun- 
damental nonlinear excitations — in the form of dark soli- 
tons. It was shown that both techniques were efficient in 
doing so. 

There are numerous interesting avenues that this work 
opens in the way of future directions. On the one hand, 
it would be useful to try to develop analytical tools to 
understand the instability of dark solitons in polariton 
condensates, as well as that of the backgrounds (i.e., the 
nodeless clouds) on which these "live". On the other 
hand, it would be interesting to extend some of the 
present considerations, such as the spectral analysis of 
stationary states, or the nonlinear interference technique 
for producing coherent structures in the 2D case to ex- 
amine some of the theoretical and potentially even ex- 
perimentally relevant (in the latter case) results thereof. 
Such studies are currently in progress and will be re- 
ported in future publications. 
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